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Abstract. We developed a time-dependent radiative model for the at- 
mosphere of HD209458b to investigate its thermal structure and chemical 
composition. Time-dependent temperature profiles were calculated, us- 
ing a uniform zonal wind modelled as a solid body rotation. We predict 
day/night temperature variations of 600K around 0.1 bar, for a 1 km/s 
wind velocity, in good agreement with the predictions by Showman & 
Guillot (2002). On the night side, the low temperature allows the sodium 
to condense. Depletion of sodium in the morning limb may explain the 
lower than expected abundance found by Charbonneau et al (2002). 



1. Introduction 

The transiting planet HD209458b should be in synchronous rotation. This im- 
plies a significant day /night temperature variation driving a strong atmospheric 
circulation with wind velocity of ~1 km/s (Guillot & Showman, 2002). Char- 
bonneau et al. (2002) constrained the sodium abundance in its atmosphere by a 
spectroscopic study. It is lower than predicted by most models which assume a 
uniform temperature field. We present a time-dependent radiative model for this 
atmosphere. It allows us, by modeling the circulation as a solid body rotation, 
to investigate the longitudinal dependence of temperature. 



2. Time-dependent radiative model 

Radiative fluxes in the atmosphere (stellar and thermal infrared) are calculated 
throught a line by line code in the spectral range 0.3-9 fim. The following opacity 
sources were included: Rayleigh scattering; H2 and He CIA absorptions; molec- 
ular species such as H2O, CO, CH4, Na, K, TiO; bound-free and H2 free- free 
absorptions. Clouds were not considered. We calculated molecular abundances 
assuming thermochemical equilibrium at each level and a solar abundance of 
the elements. Our bottom boundary condition is set by the intrinsic flux of the 
planet. Guillot & Showman (2002) showed that an equivalent intrinsic effective 
temperature of Ti Ilt =400K is required to fit HD209458b's age and radius. 
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Figure 1. Temperature at some levels as a function of longitude. 



The solution profile is calculated by solving the energy equation: 

W = ((-Zf) - (-^)) > where ~1& F */IR are respectively the heat- 

ing/cooling rates. A time-marching algorithm is used to reach a periodic solu- 
tion. We investigated the temperature as a function of longitude by modelling 
the atmospheric circulation as a solid body rotation. Fig. 1 shows the results 
for a 1 km/s wind velocity, as predicted by Showman & Guillot (2002). We 
found a 600K day/night temperature variationat the 0.1-bar level, in agreement 
with their model. A major effect is the condensation of Na on the night side. 
The sodium absorption present in the HST observations is a factor of 3 less 
than predicted by most static models (Charbonneau et al., 2002). If the sodium 
condenses in the night side, half of the limb is depleted in sodium, which might 
explain the weak HST absorption. 
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